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ABSTRACT 

Ultra-compact X-ray binaries consist of a neutron star or black hole that accretes material 
from a white dwarf-donor star. The ultra-compact nature is expressed in very short orbital 
periods of less than 1 hour. In the case of 4U 0614-1-091 oxygen-rich material from a CO 
or ONe white dwarf is flowing to the neutron star. This oxygen-rich disc can reflect X-rays 
emitted by the neutron star giving a characteristic emission spectrum. We have analyzed high- 
resolution RGS and broad band EPIC spectra of 4U 0614-1-091 obtained by the XMM-Newton 
satellite. We detect a broad emission feature at ~ 0.7 keV in both instruments, which cannot 
be explained by unusual abundances of oxygen and neon in the line of sight, as proposed 
before in the literature. We interpret this feature as O VIII Lya emission caused by reflection 
of X-rays off highly ionized oxygen, in the strong gravitational field close to the neutron star. 
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1 INTRODUCTION 

The object of our study is the ultra-compact X-ray binary (UCXB) 
4U 0614-1-091 consisting of a neutron star that accretes gas from 
a white dwarf. The sour ce is placed close to the Galactic plane 
at a distance of 3.2 kpc ( Kuulkers et al.ll2009l) towards the Galac- 
tic anticenter. The characteristic that distinguishes UCXBs from 
low-mass X-ray binaries (LMXBs) is a very short orbital period, 
it is typically less than 80 minutes. For this source, the measure - 
ments revealed an orbital period of ~ 50 min I Shahbaz et al .l2008h . 
The light curve of this object does not show signs of eclipses or 
dips, which indicates the upper limit of the inclination of the sys- 
' tem to be around 70°. The source exp eriences type I X-ray bursts 
' dSwank et alJll978l;lBrandt et al.|[l99a) and is classified as an atoll 
source iMendezetalJll997l) , which means that it shows transi- 
tions between the low/hard (island state) and high/soft state (ba- 
nana state). 

An interesting process occurring in neutron star LMXBs is reflec- 
tion. The X-ray photons emitted from a neutron star can be reflected 
by the accretion disc, leading to an X-ray emission line spectrum 
and free-free continuum. If the infalling photons are reflected by 
the innermost part of the disc, where a strong gravitational field 
is present, the observed spectral features will be broa dened by the 
relati vistic Doppler effect and gravitational redshift jpabian et al.l 
Il989l) . The emission lines can be additionally affected by Comp- 
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ton scattering in the ionized disc material, which upscatters the lin e 
photons preferentially to higher energies teallantvne et aljr2002h . 
The most prominent reflection line in the case of a disc consist- 
ing of gas with solar abundances is that of Fe Ka at ~ 6.4 keV 
and second strongest is O VIII Lya at ~ 0.65 keV. So far only 
the Fe Ka line has been reporte d in the LMXBs teeis et al.ll2009i : 
ICackett et alj2009l : lMilleil2007h . 

Evidence for reflecti on in 4U 06 14+091 was previously claimed in 
BeppoSAX data (Pir aino et alj [T999). The authors interpreted the 
spectral continuum-like feature around ~ 60 keV as the reflection 
from the optically-thick disc. A line-like feature at 0.7 keV indicat- 
ing a possible broad O VIII emission line or a set of O VII-0 VIII 
and Fe XVII-Fe XIX lines was reported in the low-resolution EX- 
OSAT an d ASCA spect ra ( Christian et al. 1994; luett et al. 200 1]). 
However. lPaerels et al.l ( 1200 ih found no evidence for em ission lines 
in the high-resolution Chan dra data. Paerels et alj j2001 ) did report 
an overabundance of neon, ljuett et alj ( 1200 ll) showed that the fea- 
ture at 0.7 keV could be eliminated in the ASCA spectra when the 
overabundance of neon is considered. However, the neutral absorp- 
tion edges of oxygen, iron and neon are unresolved in those low- 
resolution spectra. Since the 0.7 keV feature appears close to these 
edges, it is crucial for the determination of the continuum level to 
measure their depths accurately. 

We analyze simultaneous low and high-resolution spectra of 4U 
0614+091 obtained by the XMM-Newton satellite. The high- 
resolution spectra contain well resolved absorption edges, which 
allow us to improve the measurement of the abundances of oxygen. 
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Figure 1. Residuals obtained by fitting an absorbed power-law model to the M0S2 data. Top-left: Model assuming an interstellar medium (ISM) with proto- 
solar abundances of lLodde^ (2003). The region 0.6 — 0.8 k eV was excluded from the fit, Top-right: Model with an ISM where the abundance of oxygen and 
neon are allowed to vary. A similar model was proposed bv ljuett et al.l ilOOll) to explain the excess emission around 0.7 keV. The fit gives an overabundance 
of neon and an underabundance of oxygen. Bottom-left: Model with abundances of oxygen, iron and neon constrained from the RGS high-resolution spectra. 
The region 0.6 - 0.8 keV was excluded from the fit. Bottom-right: Model with parameters of the Gaussian line obtained from the RGS spectra. 



iron and neon wi th respect to the previous measurements done by 
ljuett et all ( 1200 ih . The residuals show that besides the higher than 
solar abundance of neon there is a broad emission feature present 
in the data with a peak at the position of the O VIII Lya hne. 



well as response fi les. We fit the data using the spex packagj] 
( lKaastraetal.ll996h . Errors on the fit parameters reported through- 
out this Manuscript correspond to a 68% confidence level for each 
calculated parameter (A^- = !)• 



2 OBSERVATIONS AND DATA REDUCTION 

XMM-Newton observed 4U 0614-t-091 on 2001 March 13 starting 
at 9:35 UT for ~ 11 ksec and immediately afterwards at 12:32 
UT for ~ 17 ksec. In both parts of the observation two Reflec- 
tion Grating Spectrometers (RGS) were collecting data. RGS is 
a high-resolution spectrometer covering the energy range from 7 
to 38 A (0.3 - 2.1 keV) with spectral resolution of E/AE =100 
to 500 dden Herder et al.ll200lh . In the second part of the observa- 
tion the European Photon Imaging Cameras: MOSl, M0S2 and 
pn were collecting data simultaneously with RGS. MOSl, MOS2 
and pn provide broad band spectral coverage with modest resolu- 
tion (E/AE =20 t o 50) over the energy range of 0.3 to 10 keV 
dTumer et al.l200lh . The MOS 1 and pn camera were operated in the 
Timing mode, while MOS2 was operated in the Full Frame (Imag- 
ing) mode. 

We reduce the data using the XMM-Newton Data Analysis software 
SAS version 9.0. The light curve shows no significant contamina- 
tion from soft protons. We extract source photons with pixel pattern 
equal to from the MOSl camera with rawx from 300 to 320 and 
background with rawx from 260 to 280 (beyond the source PSF). 
The net source count rate is 58.54 ± 0.06 c/s (full bandpass). We ex- 
tract the M0S2 observation with pixel pattern below or equal to 12. 
The MOS2 observation suffers from pile-up. In order to reduce the 
effect of pile-up on the source spectrum, we exclude events from 
within the circle with a radius of 24 arcsec centered on the source 
position. The background spectrum for the MOS2 observation is 
extracted from CCD-3. The net source count rate is 13.01 ± 0.03 
c/s. We extract source photons with pixel pattern less than 5 from 
the pn camera with rawx from 33 to 45 and background with rawx 
from 15 to 25. The net source count rate is 238.74 ± 0.14 c/s. 
The data collected by RGS are reduced using standard soft- 
ware pipeline which generates source and background spectra as 



3 ANALYSIS AND RESULTS 

3.1 Model with liigli neon abundance 

3.1.1 Broad band spectra 

The typical continuum model for X-ray binaries contains two com- 
ponents: a power-law, which is important in the island state and a 
black body, which is important in the banana state of the source. 
In our case adding a soft thermal component (black body) does 
not improve the fit, therefore we use only a power-law to fit the 
data. During the ob servation the source was in the island state 
dMendez et alj|2002l) . hence it is not surprising that the soft ther- 
mal component is not present. The X-ray spectrum of the source 
is affected by the i nterstellar med ium. Most of the gas in the line 
of sight is neutral dFerrierdlioOlh . therefore as to a first approxi- 
mation we use only a neutral absorption model (hot model with 
very low temperature-neutral gas limit). This model calculates the 
transmission of a plasma in coUisional ionization equilibrium. For 
the reference abundances of elements in the neutral gas we choose 
the proto-solar abundances of iLodderj ( i2003h . The data from the 
MOSl, MOS2 and pn cameras do not agree in the soft part of 
the spectrum. It was reported that the timing mode calibration 
of MOS has larger uncertainties than the imaging mode calibra- 
tion (in't Zand et al. 2008). The differences between M0S2 and 
pn spectra below 1 keV may indicate also uncertainties in the pn 
timing mode calibration. Therefore we focus on the M0S2 data 
only. In the case of the MOSl and pn data we report only the 
fit in the range 1-10 keV. The fit gives a slope of the power- 
law of 2.173 ± 0.004 with xl = 2.6 for 242 d.o.f. for MOSl and 
2.216±0.001 with;irj = 2.1 for 735 d.o.f for the pn. We fit M0S2 in 
the range 0.4- 10 keV assuming an interstellar medium (ISM) with 
proto-solar abundances. The best-fit model leaves a large residual 
in the soft X-ray part of the spectrum (Fig. 1 upper left panel). 

' http://www.sron.nl/spex 
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Figure 2. The best-fitting absorbed power-law model overlaid on the RGS spectrum extracted from the second part of the observation of the 4U 0614+091. 
The residuals show a broad emission feature with a peak at ~ 19 A. The residual at ~ 22.7 A is most probably caused by the molecular gas in the interstellar 
medium, which changes the fine structure of the absorption edge with respect to the atomic gas (see text). 



Table 1. Parameters from the best fit of the model to the MOS2 data. Pa- 
rameters fixed to the values found in the local fits to the RGS spectra are 
marked in italics. 



Parameter MOS2 





model-MOS 


model-RGS 


r 


2.23 ±0.01 


2.220 ± 0.006 




3.24 ± 0.04 


3.04 


^0 


0.50 ± 0.03 


1.0 




4.3 ± 0.2 


3.8 




1.0** 


0.23 


Eaams [keV] 




0.67 


O-Qams [keV] 




0.14 


Xl (d.o.f ) 


1.3 (279) 


1.13 (282) 



*Note the difference between the abundances obtained from fitting MOS 
and RGS data separately. 
*Parameter fixed during fitting. 

The possible solution proposed bv ljuettet all ( 1200 ll) is to allow the 
oxygen and neon abundances, whose absorption edges play a sig- 
nificant role in this region, to vary while fitting. Although the fit is 
much better, the residuals clearly show that the fit can still be im- 
proved (Fig. 1 upper right panel). In addition the fit gives an unex- 
pected significant underabundance of oxygen (see Table 1, column 
model-MOS). 



3.1.2 Constraining abundances from RGS data 

Since the spectral resolution of the M0S2 data is not sufficient to 
resolve the neon and oxygen edge, we analyse the RGS data. Since 
the source is variable, we focus on the second part of the obser- 
vation, which is simultaneous with the M0S2 observation. We use 
the slope of the power-law determined by fitting the M0S2 data as 
starting value while fitting RGS data. We do not fix the parameter 
because of the possible cross-calibration uncertainties. The high- 
resolution spectra of the source show a prominent Ne K-edge at 
14.3 A, Fe L-edge at 17.3 A and O K-edge at 23.05 A (Fig. 2). 



We measure the depth of these ab sorption edges by fit ting the spec- 
tra locally with the slab model jKaastra et id]|2002l) . This model 
calculates the transmission of a slab of gas containing individually 
selected elements with the same atomic absorption cross-sections 
as used in the neutral absorption model. The model contains the 
detailed structure of the edges. For this purpose we take the follow- 
ing energy ranges: 21 - 24 A for the O K-edge, 16 - 18 A for the 
Fe L-edge and 12 - 16 A for the Ne K-edge. For these local fits the 
power-law slope is fixed to the one found from the overall fit. The 
results are included in the Table 2. 

The abundan ce of oxygen in the ISM agrees with the proto-solar 
abundance of iLoddera l l2003 l), contrary to the results of ljuett et al.l 
f2001.) . As for the abundance of neon, we find an overabundance 
of ~ 3.8 over proto -solar, which supports the results obtained by 
IPaerels et all ( 1200 ll) . The amount of iron, that we measure is sig- 
nificantly lower than what we would expect from the ISM with 
proto-solar abundances. 



3.1.3 A feature at ~ 19 A 

We refit the full RGS spectra with the column densities of oxygen, 
iron and neon fixed to the values derived from the local fits. The 
fit gives - 2.97 for 685 d.o.f. and the residuals show structure 
around 19 A (0.7 keV). Therefore we refit the RGS spectrum ex- 
cluding the range 16-21 A. We detect a broad feature with a peak 
at 19 A in both parts of the XMM-Newton observation (Fig. 2 is 
made using only the second part of the observation). Modeling the 
feature with a Gaussian line reveals the center to be at ~ 18.6 A 
with a FWHM of ~ 4 A (see Table 2) and equivalent width of 1. 1 A 
(41 eV), which is calculated by taking the flux of the continuum 
and Gaussian line in the range 16-21 A. 

We use abundances determined from the second part of the RGS 
observation to refit the simultaneous M0S2 data with the column 
densities of each element fixed. The fit reveals an emission feature 
at around 0.7 keV (Fig. 1 lower left panel). The inclusion of the 
Gaussian line with the parameters fixed to those found in the RGS 
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Figure 3. The best-fitting absorbed power-law model with a reflection component overlaid on the RGS and MOS2 spectrum of 4U 0614+091 . The power-law 
and reflection components are additionally plotted in the dashed and dashed-dotted lines respectively. 



Table 2. Parameters from the best fit to the RGS data. 





RGS 




Parameter 


part 1 


part 2 


r 


2.37 ± 0.03 


2.26 ± 0.02 


WH[102'cm-2] 


3.2 ±0.05 


3.04 ± 0.04 


Ao 


0.93 ± 0.06 


1.01 ±0.04 




4.2 ± 0.2 


3.8 ±0.15 




0.48 ±0.13 


0.23 ± 0.07* 


xlld.o.f. 


1.82/519 


2.97/685 


Gaussian profile 


'^Gauss [A] 


18.64 ±0.17 


18.56 ±0.10 


<TGauss [A] 


4.60 ± 0.45 


3.87 ± 0.28 


EqW[k\ 


1.31 ±0.14 


1.1 ±0.1 


xlld.o.f. 


0.98/516 


1.38/682 



*Note low abundance of iron (see Sec. 4). 

data improves the fit to the M0S2 data further with;^^^ = 1.13 for 
282 d.o.f. (Table 1, Fig. 1, lower right panel). 

3.2 The broad emission feature: relativistic or Compton 
scattering broadening ? 

Whereas the Gaussian line on-top of a power-law continuum de- 
scribes the data, the model lacks a physical interpretation. A plausi- 
ble interpretation of the emission feature is that of a relativist ically 
broad ened emission line. Therefore we fit the Laor profile l lLaorl 
fl99lh to the residuals present in the RGS data. Since the peak of 
the feature is very close to the energy of the O VIII Lycy, we fix 
the wavelength of the line to 18.97 A. We fix the outer radius to 
1000 GM/c^. The fit reveals an unexpectedly high inclination of 88 
deg (Table 3). With such a high inclination, we expect to observe 
eclipses in the light curve, which are not present. It could be the 
case that the inner part of the disc, where the photons are reflected, 
has a different inclination than the binary system. The inner part of 
the disc can be warped and twisted if it is strongly irradiated by the 
neutron star in a non-uniform way (Pringle 1996). The fit gives an 
inner radius of ~ 14 GM/c^ withxl = 1-44 for 681 d.o.f., however 
we find also another solution r,„ ~ 3.5 GM/c^ with;^'^ = 1.45. 



An additional mechanism broadening the line is Compton scatter- 
ing. The Compton downscattering can be estimated using the for- 
mula cr^/ii = ET'/mpC^, where tr^ is Half Width at Half Maximum 
of the Gaussian line and E is the energy of the line in keV. The cal- 
culation gives T a 9, which seems unlikely, since it indicates an 
extremely optically-thick layer of material. This makes it impossi- 
ble that Compton downscattering alone is responsible for the line 
broadening. Therefore, the plausible solution points towards broad- 
ening by the strong gravitational field and scattering off the ionized 
material of the disc. 



3.3 Adding a reflection model 

We fit a relativistic ally broadened refle ction model r eflionx 
(Ross & Fabian 2005), available in xspec jArnaudl ll996'). which 
takes into account effects of Compton scattering in the ionized re- 
flecting material. For the relativistic broadening we use the convo- 
lution model KDBLUR (based on Laor mod el) and for the Galactic 
absorption we use the TBnew model ( Wil ms et all in prep). In or- 
der to fit the spectrum with the reflection model we need a signifi- 
cant overabundance of oxygen. In order to mimic the enhancement 
of oxygen, we decrease the abundance of iron by a factor of 5. Us- 
ing the described fit-function we find that the^- distribution around 
the (local) minimum for the inner radius does not allow for an error 
screening on the best-fit r,„ values using the current data set. Be- 
cause of additional complexity of the reflection model, we fix the 
inner disc radius to that of the innermost stable circular orbit for a 
non-rotating neutron star (6 GM/c^). The outer disc radius is again 
fixed to 1000 GM/c^. In the fit we use RGS in the same energy 
range as before and MOS2 in the range from 1.5 to 10 keV, in or- 
der to have good constraints on the power-law slope and the abun- 
dances of oxygen, iron and neon. The best fit yields a xl = 1-28 
for 1454 d.o.f. The best-fit inclination is ~ 57 deg and the best-fit 
ionization parameter is ~ 222 erg cm s"' (Table 3), which indicates 
that indeed some Compton scattering is present as well as relativis- 
tic broadening effects. 
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Table 3. Parameters from the best fit to the second part of the observa- 
tion: RGS data fitted with an absorbed power-law and a Laor profile and 
combined RGS and MOS2 data fitted with an absorbed power-law and a 
reflection model. 



parameter 


Laor profile 
RGS-part 2 


Reflection model 
RGS & MOS2 


rm 

i [deg] 

( [erg cm s" ' ] 
Xl/d.o.f. 


3.5 ±0.4 
2.24 ± 0.07 
88.0 ±0.1 

1.45/681 


14 ± 1 
2.3 ±0.1 
86.8 ±0.2 

1.44/681 


6.0* 
2.56 ± 0.03 

57.1 ±0.5 

222+22 

1.28/1454 



*Parameter fixed during the fitting. 



4 DISCUSSION 

Using high-resolution spectra we have shown that a model with 
an overabundance of neon and underabundance of oxygen cannot 
fully describe the soft X-ra y spectra of 4U 614-1-091, contrary to 
the earlier suggestions (e.g. Jjuett et alj200lh . The spectra show an 
emission line-like feature with a peak at ~ 19 A. We investigate a 
possible interpretation of this feature as relativistically broadened 
emission line of O VIII Lycf and discuss the uncertainties intro- 
duced by using only one ISM component (neutra l gas). During the 
refereeing process of this Manuscript the work of lSchulz et al.l was 
submitted, confirming the presence of this broad O VIII emission 
line using Chandra data. 

The fact that we do not observe the Fe Ka line in the broad band 
spectra indicates that the abundances in the accretion disc differ 
from solar. Optical spectra of this source show that this is indeed 
the ca se as expected from a white dwarf donor star. Nelemans e t al.l 
( l2006h find emission lines of oxygen and carbon. We detect also a 
significant overabundance of neutral neon with respect to the ISM, 
which probably originates from the source jjuett et alj|200ll) . This 
points towards a CO or an ONe white dwarf as the donor star. In 
both cases, the amount of oxygen in the disc is greatly enhanced 
with res pect to LMXBs with a main-sequence mass donor. The 
model o f iBallantvne et alj ( |2002| ) predicts that with the significant 
amount of oxygen in the disc the O VIII line becomes the most 
prominent line in the reflection spectrum. We have shown that us- 
ing a currently available reflection model it is possible to get a good 
fit by lowering the abundance of iron. The fit gives a reasonable 
value of the inclination together with significant ionization of the 
disc material. 

In our analysis we took into account only neutral gas ISM com- 
ponents. Since the source spectrum has substantial absorption of 
the Galactic ISM, we expect other components like dust, molecu- 
lar gas, warm and hot ionized gas to play a role as well. We find 
low abundance of iron, which c an be explained by t he screening ef- 
fect in the ISM dust grains (see lKaastra et alj2009l) indicating iron 
bound in dust in the line of sight. The residual at 22.7 - 23.2 A 
shows that oxygen may be bound in molecules, which significantly 
chan ges the fine structure of the edge compared to the atomic case 
jPin to et al. 2010). The narrow features close to the neon edge hint 
at the presence of ionized neon indicating warm gas, however we 
cannot constrain the amount owing to the limited signal-to-noise 
ratio of the data. The hot gas in the outer parts of the Galaxy, where 
the source is placed, originates from supernova remnants. The radio 
map of 4U 0614+091 shows no evidence of a supernova remnant 
in the line of sight (Miller- Jones, private communication). 
We attempted to constrain the inner radius of the accretion disc. 



however we obtained two solutions with comparable x~ value: 
r„, ~ 3.5 and 14 GM/c-. Deeper observations are necessary to mea- 
sure the inner radius accurately. 
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